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ABSTRACT 

Present  numerical  weather  products  do  not  include  a  current  analysis 
or  a  forecast  of  large  scale  occurrence  of  cloudiness  and  precipitation. 
The  model  used  in  this  paper  uses  vertical  velocities  and  dew-point  de- 
pressions at  the  levels  of  850,  700  and  500  mb  as  inputs  to  a  diagram 
which  determines  the  accompanying  weather  parameter.   Good  results  were 
obtained  for  a  current  analysis,  but  the  prognosis  scheme  used  proved  to 
be  unsatisfactory  at  the  present  time.   The  model  is  programmed  for  the 
Control  Data  Corporation  1604  digital  computer  to  utilize  the  operational 
output  of  machine-processed  data  and  analysis  produced  by  the  U,  S.  Navy 
Fleet  Numerical  Weather  Facility  (FNWF),  Monterey,  California. 
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1.   Introduction 

This  study  was  undertaken  to  analyse  and  forecast  areas  of  large 
scale  occurrence  of  cloudiness  and  precipitation,  including  amounts  of 
each,  using  numerical  methods.   "Large  scale  occurrence"  refers  to  the 
weather  associated  with  systems  of  the  dimensions  of  the  ordinary  cyclones 
and  anticyclones  of  the  atmosphere.   The  parameters  used  to  make  the 
analyses  and  forecasts  were  basically  vertical  motions  and  moisture 
content  of  the  atmosphere  at  the  levels  of  850,  700  and  500  mb.   Present 
vertical  motions  and  forecast  vertical  motions  were  obtained  from  U.  S. 
Navy  Fleet  Numerical  Weather  Facility  (FNWF),  Monterey,  California,   The 
present  moisture  content,  expressed  as  dew-point  depression,  was  obtained 
from  FNWF  and  moisture  forecasts  are  an  essential  part  of  the  program  used 
and  will  be  discussed  later. 

Once  the  vertical  motion  and  moisture  content  at  each  level  are 
determined,  the  results  are  averaged  and  the  Ladag  diagram,  figure  1,  is 
entered  to  determine  the  accompanying  cloudiness  or  precipitation  at  each 
grid  point.   This  diagram  is  similar  to  one  developed  by  Lewis  [1].   It 
shows  cloudiness,  showers  and  precipitation  as  a  function  of  average 
vertical  motions  and  dew-point  depression  from  850  to  500  mb„ 

The  basic  numerical  program  used  evolved  from  one  written  by  Lieutenant 
Commander  E.  M.  Carlstead,  U.  So  N„,   This  investigation  involved  checking 
and  modifying  Carlstead 's  original  program  and  then  verifying  the  results 
using  hand  analysed  cloud  and  precipitation  charts. 


2  c.   Background 

Considerable  work  has  been  done  in  the  past  to  determine  what  para- 
meters are  best  related  to  the  occurrence  of  precipitation  and  cloud 
amounts „   Smagorinsky  and  Collins  [2]  considered  the  problem  of  the 
mechanics  of  the  formation  of  precipitation  based  on  the  laws  of  fluid 
mechanics  and  used  a  computer  to  solve  these  laws  in  the  form  of  certain 
mathematical  equations.   Williams  [3]  showed  a  good  relationship  existed 
between  positive  vorticity  advection  and  precipitation,   Panofsky  [4], 
using  computed  U.  S.  Weather  Bureau  Joint  Numerical  Weather  Prediction 
(JNWP)  vertical  velocities  and  the  dew-point  depression  at  850  mb,  demon- 
strated that  the  probability  of  precipitation  increased  with  increasing 
vertical  velocities  and  decreasing  dew-point  depression.   Lewis  [1], 
making  use  of  observed  values  of  700  mb  dew-point  depression  and  JNWP 
vertical  velocities,  constructed  a  nomogram  and  divided  it  into  four 
categories:   precipitation,  cloudy,  partly  cloudy  and  clear.  Another 
approach  to  the  problem  that  showed  favorable  results  was  developed  by 
Estoque  [5].   He  used  a  graphical  technique  for  predicting  precipitation 
amounts  from  a  two-level  model  with  an  assumed  vertical  velocity  profile. 
Smebye[6]  constructed  a  model  using  the  vorticity  equation  and  the  first 
law  of  thermodynamics  to  account  for  precipitation  due  to  large-scale 
motions.   His  computed  patterns  of  precipitation  amounts  agreed  well  with 
those  observed.   In  his  model  for  predicting  precipitation  amounts  in 
Seattle  and  western  Washington,  Wilson  [7]  used  vertical  velocities, 
together  with  the  water  content  and  degree  of  saturation  of  the  air, 
Wilson's  results  gave  reasonable  estimates  of  precipitation  amounts. 
Sassman  and  Allen  [8]  observed  a  good  differentiation  between  the  occurrence 
and  nonoccurrence  of  precipitation  using  vertical  velocities  and  dew-point 


depression,  averaged  over  the  surface,  850  mb  and  700  mb  levels.   Carlstead 
[9]  has  made  numerical  forecasts  of  precipitation  and  cloudiness  using  the 
Lewis  nomogram  to  determine  the  weather  category.   His  method  shows  promise 
in  predicting  large-scale  precipitation. 


3.      The  Model 

Previous   studies   by  other  authors   have   shown  useful  correlations 
between  the   occurrence   of   precipitation  and  cloudiness  with  vertical 
motions   and  moisture   content   of   the  air.      The  vertical  velocities  were 
obtained   from  Fleet   Numerical  Weather  Facility    (FNWF),   Monterey,   Calif- 
ornia,   using  a   numerical  solution  of  a  diagnostic  co-equation  developed 
by  Haltiner,   Clarke,   and   Lawniczak    [10].      The  OD-forecasts   available  were 
in  six  hour   increments.      The  actual   numerical   program  uses  w-values  which 
are   obtained  by   the  hydrostatic   relationship 

w  =  -  —      „  (1) 

Pg 

The  moisture  content  of  the  atmosphere  is  computed  at  three  levels,  850, 
700  and  500  mb.   A  reasonable  measure  of  moisture  content  is  dew-point  de- 
pression and  this  quantity  will  be  defined  as  Tp.   The  equation  used  to 

forecast  Tp  is 

^T 

'^p(T+l)  -   ^P(T)  -  t^-^P  +  -  ^  +  Vl^^"        <2) 

On  the  right  hand  side  of  equation  (2)  the  second  term  is  geostrophic  wind 

advection,  the  third  term  vertical  advection,  and  the  fourth  term  the  pressure 

effect.   The  lapse  rate  of  dew-point  depression,  ^p,  is  assumed  to  be  dry 

adiabatic  in  nature,  and  is  approximately  equal  to  8  degrees  centigrade 

per  kilometer,  depending  on  temperature  and  pressure.  In  finite-difference 

form,  equation  (2)  then  becomes 

2        AT 
-  w  ( 
d^f 


^p(T+l)  =  ^P(T)  +  l^CfpCT)'  °(T)>  X  ■  "  <  aT  +%>!  ^'-        <3) 


In  order  to  enter  the  diagram  with  averages  of  T  and  w  the  following 
equations  are  used: 

\  '    -"3(^1  ^p850  +  \   VOO  *   =1  ■'pSOo)         <*> 


and 


w  =  .333(a2  Wg5Q  +  h^   w^^^  +  c^  ^.^^^).  (5) 


To  give  greater  weight  to  low  level  moisture  effects  the  coefficients 
used  in  equation  (4)  were  a  =  1.4j  b^  1  and  c  =  ,6.   The  tentative 
coefficients  used  in  equation  (5)  were  a^  ~  b-  ~  c„  -  1„   The  Ladag 
diagram,  figure  1,  is  then  entered  with  the  two  values  computed  from 
equations  (4)  and  (5)  and  a  value  from  0  to  600  is  obtained.   The  values 
from  0  to  600  were  selected  arbitrarily  to  assist  in  computer  computations. 
The  diagram  is  reduced  to  a  series  of  tables  and  internal  values  are  ob- 
tained by  linear  interpolation. 

If  the  value  from  the  diagram  is  a  number  greater  than  550,  rain  is 
the  forecast  weather  parameter.   To  determine  the  amount  of  precipitation 
per  hour  the  following  equation  [11]  is  used: 

P  S.3,00  (ajWgsfgSO  +  >'3''700'^700  +  '^S^OO^SOq)  '   <^> 
The  constants  a,  b  and  c  were  obtained  through  the  use  of  a  numerical 
program  called  BIMED9  [12]  using  actual  values  of  precipitation  and  the 
vertical  velocity  and  temperature  at  each  of  the  levels.   This  program 
utilizes  the  six  variables  and  using  a  best  fit  procedure  determines  the 
values  of  a,  b  and  c  in  order  to  solve  for  the  variable  P.   The  value  of 
the  constants  acquired  through  the  use  of  BIMED9  were  a„  ^  .01032,  b  ^.00097 
and  c„  =  -.00325.   Because  the  program  does  not  forecast  temperatures 
directly  the  following  hydrostatic  approximation  is  used: 

std 

From  this  approximation  the  temperature  at  each  level  becomes 

D_^--D^^.-  +  290,000 

^850  ~  ^^^'^    ^  290,000      ^   '  ^^^ 

°500-°850  +  ^11.540 
T.Qo  =  268.4  (  ^"Q  ^^Q  )   ,  (9) 

'^  411,540 

and 


°500  -  °700+2^^'^2° 
T    =  251.8  (  ^^ ^^ii^— — )   .  (10) 

^^^  256,120 

Next,  if  the  value  from  the  diagram  is  between  500  and  550,  showers 
is  the  probable  forecast  weather  parameter „   To  determine  if  showers  will 
occur  the  layer  thickness,  moisture  content  and  the  underlying  surface  is 
investigated,   Curtis  and  Panofsky  [13]  found  a  good  predictor  of  convective 
activity  by  considering  the  moisture  present  in  the  900  -  700  mb  layer  and 
the  presence  of  a  heated  surface  under  a  tropical  air  mass.   Assuming  these 
relationships  apply  elsewhere  over  the  hemisphere,  it  is  first  determined 
if  a  tropical  air  mass  is  present  by  testing  to  see  if  the  thickness  be- 
tween 1000  and  500  mb  is  5,550  meters  or  greater.   If  the  air  mass  is 
tropical  in  nature  it  is  then  determined  if  sufficient  moisture  is  present 
in  the  atmosphere  to  produce  showers.   The  criterion  for  this  is  T  must 
be  less  than  6  degrees  centigrade.   If  the  above  criteria  are  satisfied, 
it  is  then  determined  if  the  point  in  question  is  over  land  or  sea.   If 
over  sea  the  shower  equation  is  used  directly.   If  over  land  the  shower 
equation  will  be  used  only  if  the  local  time  is  between  1200  and  2000, 
The  shower  equation  [11], 

P,  =  .5(1-^),  (11) 

roughly  estimates  the  rainfall  amount  to  be  50^  of  the  precipitable  water 
in  inches  per  hour. 

Finally  the  cloud  amount  is  determined  from  the  diagram.   If  the  value 
falls  between  300  and  500  cloudy  skies  are  predicted,  if  between  100  and 
300,  partly  cloudy  skies  and  for  values  less  than  100,  clear  skies  are 
forecast. 

The  printed  northern  hemispheric  maps  available  from  this  program 

include  the  T  at  850,  700  and  500  mb  at  any  forecast  time  desired.   A 
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cloudiness  and  precipitation  value  is  also  printed  for  each  grid  point, 
This  field  is  contoured  with  the  following  scheme: 

(No  symbol)  .  .  .  .  »  .   Less  than  .3  cloudiness „ 

•  o  .  o  .  .  o  o  .  o  o  o   0.3  to  0c8  cloudiness- 
.  .  .  .  „  „   0.8  to  loO  cloudiness. 


* 
0 


Area  of  active  tropical  shower  activity  and/or 
area  of  intermittent  large-scale  precipitation, 


.  o  o  .  o  ,   Area  of  large-scale  precipitation. 
A  typical  cloudiness  and  precipitation  map  is  shown  in  figure  2,   Contours 
have  been  drawn  to  separate  the  various  weather  categories.   Figure  3  depicts 
the  geographical  area  of  the  numerical  grid  used.   It  is  noted  that  all 
large-scale  precipitation  areas  are  surrounded  by  "0"  characters.   Due  to 
an  effect  of  the  contouring  program  the  numerical  value  of  showers  must  be 
chosen  between  cloudy  and  precipitation.   This  is  rationalized  to  be  accept- 
able because  precipitation  is  often  intermittent  or  showery  at  the  edges 
of  a  large-scale  precipitation  area.   The  last  two  printed  maps  are  in 
units  of  average  precipitation  per  hour  for  the  previous  time  step  and  the 
total  precipitation  from  the  beginning  of  the  forecast  period  to  the  present 
t  ime , 

The  computations  carried  out  during  this  investigation  were  made  with 
a  Control  Data  Corporation  1604  (CDC  1604)  digital  computer,  with  a  core 
storage  capacity  of  32,768  words  of  48  bits  each.   The  input  data  were  ob- 
tained from  the  analysis  of  T  and  D  and  the  forecast  of  O)  at  the  levels 
850,  700  and  500  mb  surfaces  by  the  U.  S.  Navy  Fleet  Numerical  Weather 
Facility  (FNWF) ,  Monterey,  California.   A  simplified  flow  diagram  of  computer 
operations  is  shown  in  figure  4. 

In  order  to  utilize  available  data  and  computer  subroutines  from  FNWF, 
the  equations  were  converted  to  finite-difference  form  and  scaled  for  fixed- 
point  fractional  arithmetic. 


A   1977-point,   square-net,    octagc.ial     grid,   circumscribed  by   latitude 
ION      on  a   polar  stereographic   projection   of   the   Northern  Hemisphere  was 
employed   to   represent   the   data    fields.      The   net   spacing   is    381  km  at  60N 
where   the   projection   is   true. 

The   scaled   equations   corresponding   to  equations    (3),    (6)   and    (11) 
are  equations    (12),    (13)   and    (14)   respectively. 

A  .A  <A  A         ^cr       17   A    AT        C  7 

P  i  3600    (ajWgj^Tgjg  +  bjO^g^^^^^  +  c^^.J^^^^   2"^  (13) 


A 

T 

_i 
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P     =   -5   (1   -  T?  )  2  (14) 


The  scaling  conventions  used  in  these  equations  are 

^  17 

D  =  D  •  2'-'      (cm) 

A      9 
T  =  T   •  2     (degrees  centigrade) 


"P    P 


A 

w  =  w    • 

27 

A 
P    =    P    • 

2^1 

A 

P      =    P 
^S         ^s 

•    2 

A 

T  =  T    . 

2' 

A 

m  =  m    • 

2 

2« 
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(cm  -  sec    ) 

(inches  -  hr    ) 

(inches  -  hr  "■'■) 

(degrees  centigrade) 

(non-dimensional) 

(degrees   centigrade    -    100  m        ) 


(15) 


4.   Verification  of  Results  and  Suggestions  for  Future  Study 

Verifications  were  done  by  subjective  comparison  over  the  continental 
United  States  only.   Such  verification  of  cloudiness  is  often  difficult 
due  to  obscuration  by  fog  and  darkness.   Using  the  Lewis  nomogram  as  a 
guide,  the  lines  delineating  the  separation  of  the  different  weather 
parameters  were  shifted  in  the  Ladag  diagram  until  the  best  fit  with 
actual  conditions  resulted. 

In  the  verification  of  current  analyses j  two  precipitation  and 
cloud  depiction  maps  and  one  precipitation  amount  map  will  be  discussed. 
The  shower  areas  are  included  within  the  cloudy  category  on  the  verifi- 
cations.  Figure  5  shows  the  OOZ  5  March  1964  map<,   By  inspection  it  can 
be  seen  that  the  cloud  and  precipitation  patterns  agree  quite  well.   If 
the  shower  area  were  included  on  the  numerical  map^  the  precipitation 
areas  in  northwestern  United  States  would  be  in  much  better  agreement. 
The  actual  precipitation  map  showed  that  in  this  area  the  rainfall  amounts 
ranged  from  a  trace  to  .02  inches  during  the  previous  six  hours,  therefore 
the  maps  in  actuality  are  in  good  agreement.   Figure  6,  the  12Z  6  March 
1964  verification,  shows  a  complex  cloud-precipitation  pattern.   It  is 
felt  that  programming  a  complex  weather  pattern  such  as  this  for  the  com- 
puter would  be  difficult  indeed,  and  that  additional  parameters  besides 
vertical  velocity  and  dew-point  depression  are  perhaps  necessary.   Upon 
inspection  of  the  precipitation  area  over  western  New  York  state,  it  was 
noted  that  the  dew-point  depression  at  all  three  levels,  850,  700  and  500  mb, 
was  large,  therefore  indicating  a  partly  cloudy  sky  condition.   The  low 
level  air  trajectory  at  this  time  was  from  the  WNW  over  the  Great  Lakes 
and  it  is  felt  that  sufficient  moisture  was  picked  up  to  account  for  the 
light  snow  showers  in  this  area. 


Figure  7  shows  the  OOZ  05  March  1964  verification  of  precipitation 
amounts.   The  manual  analysis  was  done  using  total  precipitation  amounts 
in  the  6-hour  period  from  18Z  04  March  to  OOZ  05  March  1964,   The  numerical 
analysis  involves  the  computed  hourly  precipitation  rate  at  OOZ  05  March 
1964  and  multiplying  it  by  six  to  get  an  amount  for  a  6-hour  period.  As 
was  previously  mentioned  the  precipitation  amounts  in  the  northwestern 
United  States  ranged  from  a  trace  to  .02  inches  during  the  previous  six 
hour  period 3  therefore  the  verification  is  not  as  poor  as  it  first  appears 
to  be.   In  the  precipitation  area  south  of  the  Great  Lakes  j,  actual  pre- 
cipitation values  ran  as  high  as  o68  inches  during  the  previous  six  hours j 
however  the  precipitation  amount  equation  computed  a  high  value  of  only 
.38  inches.   Further  data  may  permit  the  determination  of  better  values 
for  the  constants  a„,  b  and  c_o 

Results  of  the  verifications  of  the  24-=hour  numerical  forecast  of 

cloudiness  and  precipitation  patterns j  using  the  prognosis  scheme  described 

herein,  proved  to  be  rather  poor.  Two  24=hour  prognostic  maps  are  shown  in 

figures  8  and  9.   An  initial  attempt  was  made  to  forecast  T  values  using 

only  the  geostrophic  advection  term  but  results  were  unsatisfactory  in 

regions  of  developing  high  and  low  pressure  centers.   The  forecast  o)  values 

obtained  from  FNWF  for  the  days  in  question  were  far  from  what  actually 

occurred  and  it  is  felt  that  this  is  the  largest  contributing  factor  to 

the  poor  quality  of  the  forecasts.   It  must  be  remembered  that  forecast 

w  values  are  an  essential  ingredient  to  forecast  T  values  due  both  to  the 

P 

vertical  advection  and  pressure  effect  terms  in  equation  (2).   Also  the 
forecast  w  values  are  used  directly  to  enter  the  Ladag  diagram  to  deter- 
mine the  cloudiness  and  precipitation  patterns. 

It  is  felt  that  definite  improvements  could  be  made  in  the  forecast 

model  usedo   It  was  noted  that  on  the  forecast  T  fields,  certain  areas 

P 
10 


became  either  very  wet  or  very  dry.   The  normal  five  point  Jacobian,  used 

in  the  geostrophic  advection  term  may  tend  to  over -emphasize  extreme  values 

of  T     and   certainly  can  quickly  distort   small-scale    features.      These  errors 
P 

could   be   reduced   by  using   upstream  advection  or   possibly  a   nine-point  Jacobian 

with   a   smaller   truncation  error^      Similarly,    improvement   could  also  be  made 

in   the  vertical   advection   term.      The   only  T     values   available    from  FNWF 

P 

were  the  850,  700  and  500  mb  levels.   In  computing  the  vertical  derivatives 

at  850  and  500  mb  it  was  necessary  to  use  a  forward  difference  derivative 

if  w  was  positive  at  850  mb  or  negative  at  500  mb.   In  taking  this  forward 

difference  derivative  the  assumption  is  made  that  the  T  gradient  is  constant 

P 

above  and   below  700  mb,   which  certainly  does    not   always   hold„      If  T     values 
were  available   at    1000  and   400  mb  a   better   finite-difference   approximation 
could   be   used. 
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5 ,   Conclusions 

It  is  the  authors'  opinion,  in  view  of  the  results  obtained,  that  a 
numerical  analysis  of  cloud  and  precipitation  chart,  patterned  after  this 
study  be  made  available  on  an  operational  basis «   It  is  felt  that  this 
information  would  be  of  significant  value  to  the  numerical  product  user. 

The  limitations  of  an  analysis  and  forecasting  scheme  such  as  described 
herein  are  several.   It  is  felt  that  one  of  the  larger  errors  in  the  fore- 
casts stems  from  incorrect  placement  and  magnitude  of  the  forecast  vertical 
motion  centers «   This,  in  turn,  is  a  direct  function  of  the  quality  of  the 
sea -level  and  500  mb  predictions.   S^^cond,  the   precipitation  model  depends 
upon  empirical  relationships  that  determine  precipitation  amounts,  and  these 
relationships  are  not  fully  understood.   Thirds  there  is  no  allowance  for 
the  obtaining  of  moisture  from  below  the  850  mb  level  and  in  particular 
any  increase  of  water  vapor  resulting  from  surface  evaporation.   Fourth, 
the  analysis  of  the  dew-point  depression  field  by  FNWF,  Monterey,  has 
been  shown  by  Lieutenant  Commander  D.  Howard  [14]  to  be  in  error  by  as 
much  as  6  to  8  degrees  centigrade,  even  over  dense  data  areas.   This  error 
is,  of  course,  magnified  over  sparse  data  areas.   Fifth,  the  use  of  a  one- 
hour  forecast  time  step  rather  than  the  six-hour  step  now  used  should  definitely 
improve  results.   Last,  the  effects  of  latent  heat  of  condensation  where 
precipitation  is  forming  have  not  been  included  in  any  way.   In  this  con- 
nection the  effect  on  vertical  velocity  would  perhaps  be  most  important. 

To  show  the  factor  of  machine  speed  over  the  same  job  done  manually, 
the  following  figures  are  presented:   Lewis'  technique  took  about  one-half 
hour  to  compute  the  forecast  dew-point  depression  for  one  station  at  the 
single  level  of  700  mb.   The  parcel  trajectory  was  computed  using  charts 
spaced  12  hours  apart.   The  present  machine  technique  will  produce  a  24-hour 

12 


dew-point  depression  forecast  for  1977  points  using  time  steps  of  6  hours 
at  three  seperate  levels.   This  24-hour  forecast  takes  only  10  minutes j 
a  gain  factor  of  35,586.   In  addition,  the  machine  produces  the  forecast 
in  chart  form  along  with  auxiliary  charts.   Although  improvements  are 
indeed  necessary  for  this  forecast  model,  the  time  is  near  at  hand  when 
electronic  computers  will  actually  forecast  moisture  parameters  on  an 
operational  basis. 
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Figure    1   -   Ladag  diagram, 
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Enter  diagram  with  average   vertical  motion^   w,    and  average  dew-point 
depression,    T    ,    to  determine   accompanying  weather   parameter  where 


w=   .333    (^850  +  w^oO  +  ^00^   ""^  ^p 


^  p850  p700  pSOO"^ 
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Figure  2  -  Numerical  analysis  of  cloud  and  precipitation  patterns  at 
OOZ  03  March  1964  with  manually  drawn  contourSo 
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Figure   3   -   Polar   st^reographic   projecti 
the   scale   of  figure   2„ 
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Figure   4    -  Flow  diagram  of   computer  computations 
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Figure  5  -  Manual  versus  numerical  analysis  of  cloud  and  precipitation  ' 

patterns* at  OOZ  05  March  1964. 
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Figure    6   -  Manual  versus   numerical  analysis    of   cloud  and   precipitation 
patterns   at    12 Z  06  March   1964. 
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Figure  8  -  Manual  apaiysis  ot  cionid  and  precipitation  patterns  at  OOZ  04  March  1964 

^'.■'.-rsus  ?4-b.0Jur  nume'rieal  fo'recast  from  OOZ  03  March  1964 „ 
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Figure   9   ~  Manual  analysis   o|,. cloud  and   pri--- 
'■7er3us   lA-honr  mims-tisstl  iorecas. 
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